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EXCAVATION OF THE FIRST STARS 
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ABSTRACT 

The external pollution of the first stars in the Galaxy is investigated. The first stars were born in 
clouds composed of the pristine gas without heavy elements. These stars accreted gas polluted with heavy 
elements while they still remained in the cloud. As a result, it is found that they exhibit a distribution 
with respect to the surface metallicity. We have derived the actual form of this distribution function. 
This metallicity distribution function strongly suggests that the recently discovered most metal-deficient 
star HE0107-5240 with [Fe/H]= —5.3 was born as a metal-free star and accreted gas polluted with 
heavy elements. Thus the heavy elements such as Fe in HE0107-5240 must have been supplied from 
supernovae of later generations exploding inside the cloud in which the star had been formed. The 
elemental abundance pattern on the surface of stars suffering from such an external pollution should not 
be diverse but exhibit the average pattern of numerous supernovae. Future observations for a number of 
metal-deficient stars with [Fe/H]< —5 will be able to prove or disprove this external pollution scenario. 
Other possibilities to produce a star with this metallicity are also discussed. 

Subject headings: accretion, accretion disks - stars: abundances - stars: individual (HE0107-5240) - 
supernovae: general - Galaxy: evolution - Galaxy: halo 



1. INTRODUCTION 

Though the first generation stars (so called Population 
III stars) in the Galaxy must have been formed from the 
pristine gas with no heavy elements, all the stars observed 
so far contain some amount of heavy elements such as iron 
(Fe). According to our scenario, the recent discovery of the 
most metal-deficient star HE0107-5240 with an Fe abun- 
dance as low as 1/200,000 of the solar value (Christlieb 
et al. 2002) suggests that Population III (Pop III) stars 
contained long-lived low-mass objects (Yoshii & Sabano 
1979). Thus there have survived Pop III stars even in the 
present Milky Way. To identify such stars, we have in- 
vestigated how the surface metallicity of a Pop III star 
increases by accreting gas polluted with heavy elements 
and found that HE0107-5240 is much more likely to have 
been formed from the pristine gas cloud than from a gas 
cloud that had already contained Fe (see also Thoul et al. 
2002, in the case of globular clusters and field stars). 

No stars without heavy elements have been discovered 
to date probably because stars experience the external pol- 
lution. Due to their gravity, stars can accrete the inter- 
stellar matter (ISM) polluted with heavy elements. A star 
with a mass M moving at a velocity v relative to gas ac- 
cretes gas at a rate M{v) = 1.1207r {2GM/v'^)^ pv (Bondi 
1952), where p denotes the density of the gas and G the 
gravitational constant. Thus a long-lived star moving at a 
velocity of ~ 200 km s~^ due to the gravitational poten- 
tial of the Milky Way accretes gas at a rate of the order of 
10~^^ M0yr~^. The ISM has been enriched with heavy 
elements since the first stars were formed. Hence this ac- 
cretion may increase the surface metallicity to an observ- 
able level if the surface convective layer is as shallow as in a 
metal-deficient dwarf star on the main sequence of the stel- 



lar evolution (Yoshii 1981). On the other hand, a metal- 
deficient red-giant star would reduce the external pollu- 
tion to a negligible level by diluting the accreted heavy 
elements inside its more extended surface convective re- 
gion (Fujimoto ct al. 1995; Fujimoto, Ikeda, & Iben 2000). 
The mass of the surface convective layer in a red-giant star 
becomes about two orders of magnitude larger than that 
on the main sequence. Therefore the recently discovered 
most metal-deficient red-giant star HE0107-5240 has not 
increased its surface metallicity only by accretion but has 
reached such a low value through an accretion followed by 
the dilution due to the growing of the convective envelope. 

Since the accretion rate is inversely proportional to v^, 
a much more rapid accretion can be expected when a star 
remains in its birth place, that is, in a cloud where a star 
moves with a smaller velocity relative to gas. The mini- 
mum mass that can gravitationally collapse immediately 
after the recombination of the universe indicates a typ- 
ical mass of the first clouds to be about 10^ Mq (Pee- 
bles & Dicke 1968). After contraction, a cloud formed 
a self-gravitating system in an equilibrium configuration 
with a velocity dispersion cr ^ 10 km s^^ corresponding to 
the sound speed at the recombination. Thus stars in this 
cloud move at ^/2a ~ 14 km s~^ relative to the turbulent 
gas with the same velocity dispersion and accrete gas at a 
rate ^--^8,000 times larger than those moving at a velocity of 
the order of 200 km s~^ with respect to the surrounding 
ISM. The density of the cloud being greater than the mean 
density of the Galaxy, the accretion rate is also enhanced, 
whereas the lifetime tufe (< 1 Gyr) of a cloud being shorter 
than the age of the Galaxy, the accreted mass might be 
reduced by a factor of ~20. A quantitative description of 
the external pollution by accretion will be presented in the 
next section. 
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Excavation of the first stars 
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Fig. 1. — The metallicity distribution functions (MDFs) for Pop III (solid line) and II (dashed line) stars with respect to [Fe/H]. The 
functions are normalized to unity when they are integrated over the whole range of the metallicity. The points are the MDF for Pop II dwarf 
stars (Ryan & Norris 1991). 



2. ACCRETION ONTO METAL-FREE STARS 

Prom the virial theorem, the average density p of the 
gas is given by the formula p = 3fT^/ [^-kG'^M'^). The 
accretion rate is assumed constant during the Ufetime of 

the cloud for simplicity. A star thus increases its surface 
abundance of Fc according to the formula 

[Fe/H] = log [m{v) dtX^,{t)/ (X^M^^)'^ 

-log(XFo/XH)0. (1) 

Here [Fe/H] denotes the logarithmic value of the ratio of 
Fe to hydrogen (H) relative to the solar value {^X-p^j Xyi) ^, 

X-pc{t) denotes the Fe abundance of the gas at time t, Xh 
and Mmix denote the H abundance (=0.75) and mass of 
the surface convection zone in a star, respectively. Since 
the accretion rate is the only factor that depends on the 
velocity v relative to gas, the number of Pop III stars per 
unit interval of [Fe/H] (metallicity distribution function, 
MDF) can be expressed using the distribution function 
f{v) with respect to velocity as 



dN 



d[Fe/H] 



([Fe/H]) = f{v) 



dv 



d[Fe/H] 



In 10 



vfiv). (2) 



Here stars and gas are assumed to follow the Gaussian 
distribution function with the same velocity dispersion 
(7 = 10 km s~^. Eliminating v from equations (1) and 
(2) yields the MDF. 

To describe the evolution of Fe abundance X-pc{t) of 
the gas, we use a chemical evolution model by Tsujimoto, 
Shigeyama, & Yoshii (1999) that successfully reproduces 
the distribution of halo field stars with respect to their 
Fe abundances [Fe/H] observed by Ryan & Norris (1991) 
. This model assumes that stars are formed from the gas 
compressed by the blast wave generated by a single su- 
pernova (SN) and inherit the elemental abundance pat- 
tern of the compressed gas comprising the cjecta of the 
SN. It is also assumed that the star formation was ini- 
tiated by supernova explosions of Pop III massive stars. 
The star formation ends when SN remnants overlap each 
other to heat up the cloud. Since a majority of SNe eject 



roughly 0.1 Mq of Fe per event (e.g., Hamuy 2003) and 
the mass of the ISM compressed by a single SN is of the 
order of 10^ Mq (Shigeyama & Tsujimoto 1998), it is 
predicted that most stars formed in the compressed ISM 
will have a metallicity of [Fe/H]> —4 (Tsujimoto et al. 
1999). To achieve such low values, the accretion process 
onto metal-free stars must be taken into account. This is 
consistent with the observed MDF of halo field stars be- 
longing to Pop II as shown in Figure 1 . A similar scenario 
has been proposed for the formation of globular clusters 
and Pop II field stars (Parmentier et al. 1999). As men- 
tioned above, our model requires the existence of Pop HI 
stars in a gravitationally bound cloud. The recently dis- 
covered most metal-deficient star suggests that these Pop 
HI stars include some stars with lifetimes as long as the 
age of the universe. Thus the external pollution of these 
stars must influence the current MDF of extremely metal- 
deficient stars. In this model, the mass fraction x of Pop 
HI stars is a free parameter to determine the duration of 
the star formation or the lifetime of the cloud. For a given 
mass of the cloud, this fraction must be greater than a 
certain value to form at least one massive star leading to 
SN explosion. A cloud mass of 10^ Mq requires x > 10~^ 
that leads to iufe < 700 Myr (Tsujimoto et al. 1999). Here 
the initial mass function (IMF) is assumed to be the same 
as the one derived for the solar neighbourhood (Salpeter 
1955; Tsujimoto et al. 1997). 

The mass of HE0107-5240 is estimated to be 0.8 Mq 
(Christlieb et al. 2002). Theoretical stellar evolution mod- 
els suggest that the mass range of stars that were formed 
first in the Galaxy and have now evolved to red-giant stars 
must be quite narrow around 0.8 Mq because the red-giant 
phase is brief as compared to the whole lifetime. Thus the 
MDF only for 0.8 Mq stars is considered in the follow- 
ing. Theoretical calculations (e.g., Fujimoto et al. 2000) 
show that a Pop III giant star with this mass has a surface 
convection zone with M^^ix = 0.35 Mq at present. To il- 
lustrate the maximum effect of the accretion onto 0.8 Mq 
stars, we have calculated the MDF of Pop III stars with 
this minimum fraction x (= 10~^) corresponding to the 
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longest lifetime of a cloud. The result is shown in Figure 
1 together with the MDF for later generation stars (Tsu- 
jimoto et al. 1999). These MDFs have been normalized 
to unity when they are integrated over the whole range of 
mctallicity. Here the IMF is assiimed not to vary with time 
during the halo star formation. Figure 2 shows the relation 
between v and [Fe/H] given by Equation (1). A larger frac- 
tion X of Pop III stars would lead to a shorter lifetime of 
the cloud. As a consequence, the accreted mass would be 
reduced. Thus the maximum of the MDF for Pop III stars 
would be shifted toward a smaller [Fe/H] without changing 
the shape of the distribution. Of course, the MDF for Pop 
HI stars would increase relative to that for Pop II stars 
with increasing x. Therefore Figure 1 strongly suggests 
that most of stars with [Fc/H]= -5.3 such as HE0107- 
5240 belong to Pop III. Since the shape of the predicted 
distribution is invariant to the parameter x, future obser- 
vations will be able to test the accretion model described 
above. This external pollution inevitably converges the 
elemental abundance pattern on the stellar surface to the 
uniform one represented by the average yield from SNe 
with all possible progenitor masses. The observed pattern 
in HE0107-5240 (Christhcb et al. 2002) is compatible with 
the average pattern for halo field stars, although this uni- 
formity in chemical composition is questionable (Jehin et 
al. 1999). However, the observed enhancement of the car- 
bon abundance ([C/H]= —1.3) is ascribed to a dredge- up 
of C and nitrogen to the surface taking place in an ex- 
tremely metal-poor star such as HE0107-5240 (Fujimoto 
et al. 2000; Christlieb et al. 2002) and cannot be explained 
by supernovae before the birth of this star, because a su- 
pernova supplies not only C but also oxygen, which was 
not observed. 

The short lifetime of a cloud indicates that stars were 
accreting gas while they were still unevolved dwarf stars. 
Hence a dwarf star born from the pristine gas should have 
exhibited a fairly high metallicity when the cloud was dis- 
rupted by supernova explosions at time Une- Then the star 
evolved off the main sequence to become a red-giant star 
if the stellar mass is as massive as 0.8 Mq. In the mean 
time, the surface convection zone grew and diluted the ac- 
creted heavy elements. As a result, the surface metalliciy 
has decreased by about two orders of magnitudes to be 
~ 1/100,000 of the solar value. On the other hand, less 
massive Pop HI stars still remain on the main sequence 
with a mctallicity of ~ 1/1, 000 of the solar value. 

A strong stellar wind may prevent a star from accret- 
ing the ISM. However, even if a main sequence star blew 
a wind as strong as the sun, it would hardly stop the ac- 
cretion because the typical accretion rate considered here 
is 2 X 10^^"^ M0yr~^ (Fig. 2) that is significantly larger 
than the solar mass loss rate of 2.5 x 10"^'' AfQyr"^ (Cox 
2000). Furthermore, a metal- free dwarf star with a mass of 
0.8 Mq has a surface convection zone much shallower than 
the sun leading to a less active corona in HE0107-5240 on 
the main sequence. As a result, it is likely that the wind 
from HE0107-5240 was much weaker than the sun while 
this star was accreting gas. 

3. ALTERNATIVE ORIGINS OF HEAVY ELEMENTS 

Other possibilities to produce a star with this metallicity 
will be discussed in this section. 



3.1. External pollution by a single supernova 

If a low mass Pop HI star was formed as a member of a 
binary system with a massive companion, the companion 
star might have injected a part of the ejecta into the atmo- 
sphere of the dwarf star when it exploded as a SN. Since 
the average speed of ejecta is ~ 3,000 km s""'^, the mass 
injection into the dwarf star is expected to be extremely 
inefficient. If the ejecta that expand into the geometrical 
cross section of the star with the radius R are injected, 
the total mass Minj of the injected ejecta will be given 

by the formula Minj = M^j {R/2D)^, where Mej denotes 
the total mass of the ejecta and D the separation of the 
stars. This will give the maximum injected mass because 
the ejecta move at a velocity much greater than the es- 
cape velocity of the star. To explain the observed value 
of [Fe/H]= —5.3 on the surface of a star with the solar 
radius, the separation at the time of SN explosion must be 
closer than ~ 10-^^ cm if the SN ejected ~ 0.1 M© of Fe. 
In reality, the ejecta with sufficiently high velocities can- 
not be injected into the star even if they initially proceed 
toward the star. Since the radius of a red super- giant star 
could exceed this separation, these two stars would merge 
into a single massive star before SN explosion. 

If the companion star was so massive ( ^ 100 Mq) that 
it underwent a so-called pair instability supernova that 
ejected more than 10 Mq Fe (e.g., Heger & Woosley 2002), 
it was able to inject enough Fe into the dwarf star to ex- 
plain the observed metallicity even if the separation was 
as large as ~ lO^"' cm. According to c;m-rent theoreti- 
cal models for this type of supernovae (Heger & Woosley 
2002), however, it seems that no SN can inject heavy el- 
ements with the abundance pattern observed in HE0107- 
5240. Since different elements expand at different veloci- 
ties in SN ejecta, detailed hydrodynamic calculations are 
needed to explore the abundance pattern of heavy elements 
injected by this process. A pair instability supernova is 
so energetic that it can disrupt a cloud with a mass of 
10^ Mq. As a result, other Pop HI stars inside the same 
cloud would not accrete the gas. Pair instability super- 
novae might prevent most of metal-free stars from an ex- 
ternal pollution. 

3.2. HE0107-5240 as a Pop II star 

There might still be a possibility that stars with [Fe/H] = 
—5.3 were formed from gas into which the precursor SNe 
had already injected heavy elements, that is, belong to Pop 
II. It has been suggested that the elemental abundance 
pattern of a metal-deficient star (with [Fe/H] ^ — 2) was 
imprinted by a single SN and coincides with that in the 
gas eventually swept up by the blast wave generated from 
this SN (Shigeyama & Tsujimoto 1998). Thus the amount 
of Fe corresponding to [Fe/H]= —5.3 suggests that the 
precursor SN ejected only ~ 0.0005 Mq of Fe. In fact, 
the current SN observations have revealed that each of 
the least luminous SNe ever observed, SN 1997D and SN 
1999br, ejected a small amount of Fe, while a majority of 
SNe eject an amount of Fe per event at least one order of 
magnitude larger (see Hamuy 2003; Zampieri et al. 2003). 
Theoretical SN models suggest that the deficiency of Fe 
must be a result of the fall back of some Fe in the ejecta 
onto the compact remnant at the center (Turatto et al. 
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Fig. 2. — The mass accretion rate (solid line; in the unit of Mq yr ^) and the velocity of a star relative to gas (dashed line) as functions 
of the surface metallicity [Fe/H] of Pop III star with a surface convection zone of mass Mmix = 0.35 Mq. 



1998). As a consequence, the ejecta exhibit an apparent 
enhancement of a-elements such as magnesium (Mg) rela- 
tive to Fe. By contrast, the observed ratio in HE0107-5240, 
[Mg/Fe]= 0.2, is comparable with the average value for 
halo field stars. Therefore it is unlikely that HE0107-5240 
inherited heavy elements ejected by a subhiminous SN. On 
the contrary, two extremely metal-deficient stars CS22949- 
037 and CS29498-043 exhibit distinct enhancements of a- 
elements with respect to Fe such as [Mg/Fe]> 1 (Depagne 
et al. 2002; Aoki et al. 2002). These stars might have in- 
herited heavy elements from subluminous SNe (Tsujimoto 
& Shigeyama 2003). 

There is, of course, inhomogeneity in the matter swept 
up by a single SN. Thus it could be expected that stars 
with [Fe/H]=-5.3 were formed even from luminous SNe. 
A numerical study of metal enrichment by the first SN in 
an inhomogeneous ISM apparently denies this possibility. 
3-D hydrodynamic calculations of a SN explosion at differ- 
ent sites has shown that there is little gas with [Fe/H]=- 
5.3 inside the matter swept up by the SN (Nakasato & 
Shigeyama 2000). In addition, the MDF must be easily 
distinguished from that presented in Figure 1 because the 
mechanism to acquire Fe are quite different from the accre- 
tion. The MDF for stars with [Fe/II]< —5 is indispensable 
for identifying the origin of heavy elements in these stars. 
Further differences must exist in the abundance pattern 



of each star that reflects the precursor SN and exhibit a 
feature different from each other. This is in contrast with 
the consequence of the accretion mentioned above. 

4. CONCLUSIONS 

Here we have discussed the gas accretion onto metal- 
free stars while they stayed in a cloud with a mass of 
~ 10'' Mq and predicted the MDF for Pop III stars as 
a consequence of the external pollution. After the disrup- 
tion of the cloud, stars have been moving at ~ 200 km s~^ 
under the gravitational potential of the whole galaxy and 
no longer accreting a substantial amount of matter. This 
MDF strongly suggests that the recently discovered most 
metal-deficient star HE0107-5240 belongs to Pop III. The 
expected elemental abundance patterns on the surfaces of 
these Pop III stars with external pollution should not be 
diverse but exhibit the average pattern of numerous SN ex- 
plosions. The observed pattern in HE0107-5240 satisfies 
this requirement. The argument presented here can be 
tested by future observations for a number of extremely 
metal-deficient stars with [Fe/H]< —5. 
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